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Abstract 



Cluster production is investigated in central collisions of Ca+Ca, Ni+Ni, Zr+ Zr, 
96 Ru+ 96 Ru, Xe+Csl and Au+Au reactions at 0.4A GeV incident energy. We find 
that the multiplicity of clusters with charge Z > 3 grows quadratically with the 
system's total charge and is associated with a mid-rapidity source with increasing 
transverse velocity fluctuations. When reduced to the same number of available 
charges, an increase of cluster production by about a factor of 5.5 is observed in 
the mid-rapidity region between the lightest system (Ca+Ca) and the heaviest one 
(Au+Au). The results, as well as simulations using Quantum Molecular Dynamics, 
suggest a collision process where droplets, i.e. nucleon clusters, are created in an 
expanding, gradually cooling, nucleon gas. Within this picture, expansion dynamics, 
collective radial flow and cluster formation are closely linked as a result of the 
combined action of nucleon-nucleon scatterings and the mean fields. 



In energetic central heavy ion collisions it is generally assumed that, after go- 
ing through an early stage of hot and compressed nuclear matter, the system 
undergoes a substantial expansion causing local cooling before freezing out. 
At beam energies above 10A GeV the current picture is that the hot sys- 
tem while cooling passes from a phase involving at least partially deconfined 
quarks and gluons into a purely confined hadronic phase. At energies below 
1A GeV the hot phase is still predominantly a nucleonic gas which, however, 
in the expansion phase can partially 'liquefy' i.e. clusterize in analogy to the 
processes used in clusterization devices for atomic physics [1], although under 
less controlled conditions. In both energy regimes the aim is to determine ba- 
sic parameters, such as the critical temperature T c or the latent heat of the 
(first order) phase transition. Due to the finite size of the nuclear systems 
available in accelerator physics and the complexity of the dynamics of heavy 
ion reactions, convincing progress on this frontier has proven to be a diffi- 
cult task. Concerning the liquid-to-gas transition onsets of plateau's in caloric 
curves [2] have been interpreted [3] in ways that relate indirectly to first order 
transitions. More direct signatures [4,5], such as negative heat capacities [6] 
have been subjected to critical reviewing [7]. Claims to the determination of 
T c [8,9,10] vary in the proposed values and require model assumptions, that 
seem to be in conflict with some experimental data [11]. Indications for the 
existence of a negative compressibility region in the nuclear phase diagram, 
leading to spinodal instabilities characterized by enhancement of events with 
nearly equal-sized fragments, have been obtained [12]. The very small cross 
sections for this phenomenon were justified with microscopic simulations. One 
of the key assumptions in many works is that multifragmentation is a unique 
mechanism related directly to subcritical and/or critical phenomena [5]. In 
recent theoretical simulations using Nuclear Molecular Dynamics (NMD), to- 
gether with a backtracing method, Bondorf et al. [13] have argued that two 
mechanisms coexist: a dynamical rupture of spectator-type, relatively cold, 
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fragments and a second process where nucleon-nucleon collisions generate the 
seeds for completely new fragment creation, coalescing nucleons which were 
initially far from each other in phase space. 

In this Letter we present data for central collisions in symmetric heavy ion 
systems at 0.4A GeV that strongly support these theoretical findings [13]. In 
order to better understand the finite-size problem, we have varied the sys- 
tem size from Ca+Ca (Z sys = 40) to Au+Au (Z sys = 158), investigating five 
systems of different size. We find that the multiplicity of heavy clusters with 
charge Z > 3, when reduced to the same number of available charges, grows 
linearly with system size and is associated with a mid-rapidity source with 
increasing transverse momentum fluctuations. This is in strong contrast to 
multifragmentation of quasi-projectiles in the Fermi energy regime or of spec- 
tator matter at SIS/BEVALAC energies (0.1-2A GeV) where a high degree 
of 'universality' [14,15] was observed, with among others an apparent system- 
size independence. As parallel studies of the same systems [16] have shown 
an increasing degree of stopping as well as increasing generation of flow when 
passing from light to heavy systems, we can conclude that the increased flow 
leads to a cooling process favouring gradual 'liquefaction'. This is a non-trivial 
finding: from earlier experimental studies [17] a suppression of heavy fragment 
production in systems with strong collective expansion was inferred and the- 
oretical works [18] have predicted that strong flow gradients would prevent 
coagulation. A similar behavior is seen in two-particle correlations where the 
effective radii of homogeneity are diminished by flow [19]. 

The data were taken at the SIS accelerator of GSI-Darmstadt using various 
heavy ion beams of OAA GeV and the large acceptance FOPI detector [20,21]. 
In the experiments involving the systems 40 Ca+ 40 Ca, 96 Ru+ 96 Ru, 96 Zr+ 96 Zr, 
and 197 Au+ 197 Au, particle tracking and energy loss determinations were done 
using two drift chambers, the CDC (covering laboratory polar angles between 
35° and 135°) and the Helitron (9° — 26°), both located inside a superconduct- 
ing solenoid operated at a magnetic field of 0.6T. A set of scintillator arrays, 
Plastic Wall (7-30)°, Zero Degree Detector (1.4°-7°), and Barrel (42°-120°), 
allowed to measure the time of flight and, below 30°, also the energy loss. All 
subdetector systems have full azimuthal coverage. Use of CDC and Helitron 
allowed the identification of pions, as well as good isotope separation for hydro- 
gen and helium clusters in a large part of momentum space. The identification 
of heavier clusters (Z > 3), by nuclear charge only, was restricted to the polar 
angles covered by the Plastic Wall and the Zero degree detector. In a second 
setup the Helitron was replaced by an array of gas ionization chambers [20], 
the PARABOLA, allowing charge identification of heavier clusters up to nu- 
clear charge Z = 12. The systems 58 Ni+ 58 Ni, 129 Xe+CsI, 197 Au+ 197 Au were 
studied in this experiment. The data for Au on Au, measured with both se- 
tups, were found to be in excellent agreement. Further details on the detector 
resolution and performance can be found in [20,21]. 
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Collision centrality selection was obtained by binning distributions of the ratio, 
Erat [22], of total transverse and longitudinal kinetic energies. In terms of the 
scaled impact parameter, b^ ' = b/b max , we choose the same centrality for all 
the systems: b^ < 0.15. We take b max = 1.15(Ap + Aj. 3 ) as effective sharp 
radius and estimate b from the measured differential cross sections for the 
Erat distribution using a geometrical sharp-cut approximation. In this energy 
regime Erat selections show better impact parameter resolution for the most 
central collisions than multiplicity selections [22,23] and do not imply a priori a 
chemical bias. Autocorrelations in high transverse momentum population, that 
are caused by the selection of high Erat values, are avoided by not including 
the particle of interest in the selection criterion. 

The analysis of the particle spectra involves some interpolations and extrap- 
olations (30% in the worst case) to fill the gaps in the measured data. The 
two-dimensional method used to achieve this has been extensively tested with 
theoretical data generated with the transport code IQMD [28] applying appa- 
ratus filters. Details will be published elsewhere [29]. 

Since the present study is restricted to central collisions of symmetric systems, 
we require reflection symmetry in the center of momentum (c.o.m.) and, use 
azimuthally averaged data. Choosing the c.o.m. as reference frame and orient- 
ing the z-axis in the beam direction the two remaining dimensions are charac- 
terized by the longitudinal rapidity y = y z , given by exp{2y) = (1+/3 Z )/ (1—/3 Z ) 
and the transverse component of the four-velocity u, given by u t = j3 t ^f. Fol- 
lowing common notation /3 is the velocity in units of the light velocity and 
7 = 1/J\l — (3 2 ). Later we shall also use the transverse rapidity, y x , which is 
defined by replacing f3 z by (3 X in the expression for the longitudinal rapidity. 
The x-axis is laboratory fixed and hence randomly oriented relative to the re- 
action plane, i.e. we average over deviations from axial symmetry. Throughout 
we use scaled units y^ = y/y p and = u t /u p , with u p = (3 p "y p , the index p 
referring to the incident projectile in the c.o.m.. 

An example of a reconstructed distribution for emitted Li ions in central col- 
lisions of Au on Au is shown in Fig. 1. The lower part of this figure illustrates 
the result of the extrapolation to An using the two-dimensional fit method. 

Turning now to the presentation of results, we show in Fig. 2 charged particle 
multiplicity distributions as a function of nuclear charge Z. The data for four 
out of the six measured systems are plotted. As the surface of nuclei has a 
finite thickness one expects some degree of transparency even in collisions 
with perfect geometrical overlap. To minimize such 'corona' effects, we show 
midrapidity data, \y^\ < 0.5. As observed earlier [22], we find that the heavy 
cluster [Z > 2) multiplicities, M hc , decrease exponentially with the nuclear 
charge, i.e. M hc (Z) ~ exp(— c* hc Z), however, the slope parameter c* hc is seen to 
vary with the system-size (see Fig. 2 and Table 1): heavy cluster production is 



4 



Au+Au 400A MeV 

i | i | i | i | i | r 




-1.0 -0.5 0.0 0.5 1.0 



I I I I I 


i i i i i 

Z=3 







L_ 

-1.0 -0.5 0.0 0.5 1.0 



y (0) 

Fig. 1. Invariant distribution dN / {uf^ du^ dy^>) of Li ions emitted in central col- 
lisions of Au+Au at 0AA GeV. The various grey (color) tones correspond to yields 
differering by a factor of 1.5. The upper panel shows the measured data includ- 
ing the use of reflection symmetry and interpolations with the two-dimensional fit 
method. The lower panel shows the extension to 4ir. 

favoured for larger systems. 

In the present context we refrain from calling heavy clusters (he), Z > 2, 
'intermediate mass fragments' (IMF), because in this energy regime they are, 
actually, the heaviest fragments ('droplets') in central collisions. Charge bal- 
ances show unambiguously that there is no heavy remnant [Z > l/6Z sys ) with 
a sizeable (> 1%) probability (that one might be tempted to call 'liquid'). The 
slope parameters, obtained for the range Z = 3 — 6, are listed in Table 1 both 
for the mid-rapidity data, c* hc , and the Air data, Ch c - For Ca+Ca reliable data 
beyond Z = 4 could not be obtained. The value Ch c = 1.224 ± 0.043 for 
Au+Au can be compared with our earlier [22] value of 1.170 ± 0.018 which 
was obtained from a fit over a larger Z range. Although a different method 
to extrapolate to 47r was used, the main reason for the modest deviation can 
be traced to the somewhat higher centrality achieved with the present setup 
compared with the 'PHASE I' setup used in our earlier work, which covered 
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Fig. 2. Charged particle multiplicity distributions - dN/dZ - measured in central 
collisions (b^ < 0.15) of Au+Au, Xe+Csl, Zr+Zr and Ni+Ni reactions at an inci- 
dent energy of OAA GeV. The data are presented for c.o.m. rapidities : \y^\ < 0.5. 
The straight lines are exponential fits to the data in the interval 3 < Z < 6. Errors, 
if not shown, are smaller than the symbol sizes. 

only polar angles below 30°. 

If one were to interprete the slopes as being an indicator of a global freeze-out 
temperature of systems in chemical equilibrium, the qualitative conclusion 
would be that smaller systems appear to be hotter. We recall that all sys- 
tems are studied at the same incident energy and for the same centrality, 
= b/b max < 0.15. However, in complex reactions leading to many outgo- 
ing channels with an apparently simple (exponential) statistical distribution, 
integrated charge distributions give limited information and therefore do not 
allow to draw convincing conclusions on the possible emergence of a final state 
in equilibrium. 

More details of the mechanism at work are seen in Fig. 3 where we com- 
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Fig. 3. Scaled longitudinal (left) and transverse (right) rapidity distributions of 
single protons (top) and Li clusters (bottom) measured in central collisions for 
different symmetric reactions. The various system charges Z sys are indicated in the 
Figure. The ordinates are normalized to a common reference system with charge 

Z sys = 100. 

pare the evolution of the longitudinal and the transverse rapidity distribu- 
tions (left and right panels, respectively) with system size for the 'gas' (pro- 
tons) which is prevalent for the light (surface-dominated) systems, and the 
'droplets' which appear in increasing numbers at mid-rapidity as the system- 
size increases (lower left panel) , accompanied by a broadening of the transverse 
velocities (lower right panel). We plot the smoothened data emphasizing the 
gradual evolution with system size. The statistical errors are smaller than the 
symbol sizes, while the absolute systematic errors are 10%. Relative errors are 
expected to be smaller by a factor of two or more. 

Note that all yields in Fig. 3 are reduced to a constant system size of 100 
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Table 1 

Characteristics of heavy clusters production in central heavy ion collisions at OAA 
GeV. 



System 


Che 


r* 
tie 


M hc 


M he 


vart(Li) 








±10% 


±10% 


±0.02 


40 Ca + 40 Ca 






0.42 


0.133 


0.094 


56 Ni + 40 Ni 


1.44 ± 0.07 


1.63 ± 0.08 


0.75 


0.32 


0.117 


96 Zr + 96 Zr 


1.441 ± 0.046 


1.49 ± 0.06 


1.55 


0.84 


0.139 


96 Ru + 96 Ru 


1.415 ± 0.042 


1.37 ± 0.06 


1.65 


0.89 


0.144 


108 Xe + 133 Cs 12 7l 


1.276 ± 0.040 


1.336 ± 0.044 


2.29 


1.35 


0.147 


197 Au + 197 Au 


1.224 ± 0.043 


1.217 ± 0.044 


4.66 


2.99 


0.177 



nuclear charges by multiplying with (100/Z sys ), where Z sys is the total system 
charge. Closer inspection of the system-size evolutions in these reduced scales 
reveals two other noteworthy features: 1) The 'missing' proton vapour in the 
heaviest system, that served as source of the developping clusters, is limited to 
smaller transverse rapidities, \y£'\ < 0.7, see the upper right panel and note 
the enlarged abscissa scale in the lower right panel. 2) Except for the Ca+Ca 
system, heavy cluster production is system-size independent around longitu- 
dinal rapidities, \y°\ > 0.8, see the lower left panel. This latter point is not 
trivial, although 'universal' features in the partition of excited spectator mat- 
ter are well established [14]. In general, however, the collisions investigated in 
ref. [14] were more peripheral and 'spectators' could be clearly identified as rel- 
atively narrow peaks in the longitudinal rapidity distributions near projectile 
(or target) rapidity. In the present cases of high centrality, spectator matter, 
if any, is not readily isolated in the measured rapidity distribution, except for 
the two-peak structures in the Ca±Ca reaction. But even in that case the 
evolution towards mid-rapidity is continuous. When we select more peripheral 
collisions 'spectator' peaks appear in all systems and we observe [24,25] at this 
energy the same universal features as in ref. [14]. 

A way to summarize the observations and at the same time confronting them 
with theoretical simulations is presented in Fig. 4. Starting with the lower left 
panel, we show the system-size dependence of the rapidity-integrated multi- 
plicity, M hc , of heavy clusters. Although we plot again reduced yields (to 100 
incoming protons), we observe a remarkable linear increase: the least squares 
fitted straight line follows the data with an accuracy of 4% (as mentioned 
earlier, the relative accuracy of the data points is expected to be better than 
the indicated 10% systematic error bars). This means that the system-size 
dependence of M hc has a term proportional to Z 2 sys or A 2 sys (the system mass 
squared) . 
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There is a small irregularity associated with the two data points due to Zr+Zr 
(Z sys = 80), resp. Ru+Ru (Z sys = 88) which represents systems with the same 
total mass (A sys = 192), suggesting that besides the size dependence there is 
also an isospin dependence of heavy cluster production. The effect is however 
at the limit of the error margins and therefore will not be discussed any further 
in the present work. 

Besides the A 2 sys term in Mh c (A sys ) , a second piece of information comes from 
looking at multiplicities, M£ c , confined to the midrapidity interval (\y^\ < 
0.5) (see lower left panel of Fig. 4): the restricted data run parallel to the in- 
tegrated data, showing that the quadratic term is associated virtually entirely 
with the midrapidity region. For comparison we also show the same midrapid- 
ity data before acceptance corrections. While the corrections are important, 
and enhance the effect, it is clear that they are not creating the effect. In the 
midrapidity region the reduced cluster production is about a factor 5.5 higher 
for Au+Au than for Ca+Ca. We note that for the lightest system, Z sys = 40 
is still large compared to the most abundant [Z = 3) heavy cluster. Of course 
the linear trends cannot go on indefinitely as the number of heavy clusters per 
hundred protons cannot exceed 33 by definition. We are far from this trivial 
limit, however. 

We have also reduced M£ c by dividing by the sum of charges accumulated at 
midrapidity (\y^ < 0.5), Z mi d y , rather than by the total system charge, Z sys . 
In this case the straight line fitted to the data (not plotted to avoid overloading 
the Figure) passes very close to the origin of the axes: the relationship M£ c x 
100/Z midy = 0.0342Z midy reproduces the six data points with an accuracy of 
6%. This means that, well within our systematic errors, we can say that all 
the mid-rapidity heavy cluster production rises with the square of the number 
of 'participant' nucleon number. 

The associated observation of increasingly broad transverse rapidity distribu- 
tions as the system size is increased, shown in Fig. 3, is summarized in the 
upper left panel of Fig. 4 where the variances var(y^') of these distributions 
(taken for data within \y^\ < 1) are plotted and also seen to rise linearly 
with the system size, although with an offset. The observables M^ c , M£ c and 
var(y^) are also listed in Table 1. In a global thermal equilibrium picture the 
variance of rapidity (velocity) distributions is a measure of the temperature 
and hence, in this scenario, one would conclude that the heaviest system is 
the hottest system, in contradiction to our earlier conjecture, inferred from the 
observed partitions, that the lightest system is the hottest. A way out of this 
contradiction is to introduce radial flow: generated during the expansion, it is 
accompanied by a local cooling. This mechanism is seen to be stronger for big- 
ger systems. Radial flow, assessed by modelling hydrodynamic expansion [26], 
or treated phenomenologically [22] , has been introduced by our Collaboration 
with some success in describing the deviations of momentum space distribu- 
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Fig. 4. Summary of heavy cluster data measured with FOPI (left panels) and cal- 
culated using IQMD (right panels). Lower panels: Average multiplicities of heavy 
clusters per 100 incoming protons (i.e. Mh c x W0/Z sys ) as a function of the system 
charge Z sys . In decreasing order the data represent rapidity- integrated (Air) data 
(open triangles), data confined to the midrapidity interval \y^\ < 0.5 (full dia- 
monds) and, in the FOPI case, data limited to the acceptance of the apparatus, i.e 
before acceptance corrections (open circles). Upper panels: Variance of the trans- 
verse rapidity distributions for Li ions versus system charge. All straight lines are 
linear least squares fits to the respective data. Note that the IQMD multiplicities 
are multiplied by a factor three. The meaning of the marked square symbols in the 
IQMD data is explained in the text. 

tions from a global thermal scenario, notably the fragment mass dependences. 
But some inconsistency in accounting for the cluster yields [22] remained. 
This difficulty could be due to the failure of correctly accounting for non- 
equilibrium effects caused by the coronae of the nuclei and, more generally, by 
partial transparency [16]. 

In principle, such non-equilibrium effects can be handled by simulation codes 
based on transport theory [27]. We have used the code IQMD [28] based on 
Quantum Molecular Dynamics [30], QMD, to see if we could reproduce the 
features of our data. Clusters are identified after a reaction time of 200 fm/c 
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using the minimum spanning tree algorithm in configuration space with a co- 
alescence radius of 3 fm. For each system-energy we have generated 50000 
IQMD events over the complete impact parameter range which were subse- 
quently sorted using the Erat criterion to obtain an event class with compa- 
rable centrality to that selected in the experiment. The reduced heavy cluster 
multiplicities (4tt integrated or at mid-rapidity) versus system size are plot- 
ted in the lower right panel of Fig. 4. As in our earlier studies [22], we find 
that IQMD, as well as other realizations of QMD [31], underestimates clus- 
ter production: the yields in the Figure are multiplied by a factor of three. 
However, qualitatively, the experimental size dependence is reproduced, in- 
cluding the dominant effect of mid-rapidity emissions in accounting for the 
linear rise. The associated rise of the variance of the transverse rapidity with 
system size is reproduced almost quantitatively. In this context it is useful 
to realize that the nuclear stopping phenomenon, as quantized by the shape 
of proton and deuteron longitudinal rapidity distributions is rather well de- 
scribed by IQMD [32] at incident energies around 400A MeV. For heavy clus- 
ters the simulations also reproduce the yields near projectile (target) rapidity: 
the underestimation concerns mid-rapidity emissions and might be a general 
deficiency of semiclassical approaches such as IQMD that tend to converge 
towards Boltzmann statistics after many single nucleon collisions (leading to 
mid-rapidity population) rather than conserving the initial Fermi-destribution 
of nucleons. 

To shed some light into the cluster creation mechanism, we have performed, 
for Au+Au, a calculation where the elementary nucleon-nucleon cross sections 
were raised by a factor two. The results (full squares marked '2<r' in the right 
panels) show an increased transverse rapidity variance, as one would expect 
from the increase of the elementary collision frequency, but also a rise of the 
heavy cluster multiplicity, shown for the interval ly^l < 0.5, supporting our 
observation that cluster production is in a direct correlation with flow. The 
introduction of more copious binary interactions increases the adiabaticity 
(and hence the cooling effect) of the expansion in a twofold way: 1) the local 
equilibration is faster and 2) the expansion is slower (since the diffusion time 
is larger). Clearly, we have found evidence for the 'second process' suggested 
in [13] where nucleon-nucleon interactions act as seeds for self-organization 
leading to new clusters. 

A key question is whether cluster production is sensitive to relevant features 
in the nuclear phase diagram, in particular to the existence and location of a 
critical point and of liquid-vapor coexistence curves. When using the IQMD 
code with its two available momentum dependent Equations of State, EOS, as 
input, a so-called 'stiff' one with an incompressibility K = 380 MeV around 
saturation density, or alternatively a 'soft' EOS, K = 200 MeV, we find sen- 
sitivity of cluster production. While all calculations presented so far in Fig. 4 
were done with the soft EOS, one calculation, again for Au+Au, was performed 
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assuming a stiff EOS: see the open squares (marked 'HM') in the right panels 
of Fig. 4. The variance var(y^') is little changed, i.e. it is caused primarily 
by the nucleon-nucleon collision frequency which is only indirectly affected 
by switching to another (mean field) EOS. Cluster production in the interval 
\y^\ < 0.5, however, is lowered when assuming a stiff EOS. 

It is too early to draw specific conclusions from this potentially interesting find- 
ing. What is missing, so far, is the underlying phase diagram implied by IQMD 
in its various options. Recently it was shown [33] with codes based on Anti- 
symmetrized Molecular Dynamics (AMD) that caloric curves can be predicted 
in controlled scenarios (fixed volume or pressure). This presents an important 
and necessary link between microscopic models for equilibrium thermodynam- 
ics (and nuclear structure as well) and the complex dynamic situation found 
in heavy ion collisions. 

In summary, cluster production has been investigated in central collisions 
of Ca+Ca, Ni+Ni, 96 Zr+ 96 Zr, 96 Ru+ 96 Ru, Xe+Csl and Au+Au reactions at 
OAA GeV incident energy. We find that the multiplicity of clusters with charge 
Z > 3, when reduced to the same number of available charges, grows linearly 
with system size and is associated with a mid-rapidity source with increasing 
transverse velocity fluctuations. An increase by about a factor of 5.5 is ob- 
served in the mid-rapidity region between the lightest system (Ca+Ca) and 
the heaviest one (Au+Au). The results, as well as simulations using Quantum 
Molecular Dynamics, suggest a collision process where droplets are created in 
an expanding, gradually cooling, nucleon gas. Expansion dynamics, collective 
radial flow and cluster formation are closely linked resulting from the combined 
action of nucleon-nucleon scatterings and the mean fields. 

Finally, we note that global stopping and directed sideflow data are available 
for the same systems [16]. The simultaneous reproduction of both 'repulsive' 
observables (sideflow) and 'attractive' observables (radial flow and its effect on 
the degree of clusterization) with the same transport code will be a challenging 
task. The reward, hopefuly, will be a more precise mapping of the nuclear phase 
diagram, including the liquid- vapour transition, than was possible in the past. 

This work was supported in part by the French-German agreement between 
GSI and IN2P3/CEA (project No 97-27). 
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